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Recurrent swelling of horizontally shaken granular material 
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This Letter reports an experimental study of the dynamics of granular materials when vibrated 
horizontally. Convective motion is observed over wide range of parameters, similarly to the known 
effect observed under vertical vibration. However, under horizontal vibration we observe a striking 
novel effect on the oscillating height of the surface of the material: the frequency of the surface 
oscillations are totally decoupled from the frequency of the drive. We explain the effect as resulting 
from an interplay between Reynolds dilatancy due to convective motion and mechanical stability of 
the material. 

PACS: 46.10 +z, 83.70 Fn, 24.60.-k. 
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When a rectangular container with granular material 
is shaken either vertically or horizontally one observes 
convection. For the case of vertical shaking this effect 
has been intensively studied experimentally, analytically 
and by means of computer simulations (s. e.g. ||l|,g,^). 
Many results on vertical shaking are known and there 
is a certain degree of understanding of these phenom- 
ena. For horizontal shaking convection cells have been 
reported only recently Mpl and presently one is far from 
understanding these complex phenomena. 

In the present letter we report a new effect experimen- 
tally observed in horizontally shaken granular material. 
As far as we know this is the first dynamical effect where 
the frequency of forcing decouples from the frequency 
of the effect, i.e. the system reveals an inherent time 
scale Q. 

The new effect, which we want to focus here, is re- 
current swelling of horizontally shaken material: For a 
certain range of amplitude A and frequency / of the si- 
nusoidal motion of the container x{t) — A cos(27r / 1) and 
for a certain region of filling level we observe a time de- 
pendent variation of the volume of the material. We will 
show that the fluctuations of the height of the material, 
the intensity of the convection flow and the energy dis- 
sipation rate in the flowing material are closely related 
effects. 

There are several experimental investigations of ex- 
citing time-dependent phenomena in vertically vibrated 
containers Q,^- While the time scales of these effects 
are comparable with the period of forcing, for horizontal 
vibration we find a main difference: The effect reported 
in this letter has a period which is up to several hundred 
times larger than the period of oscillation f~^ . Hence we 
do not observe a period doubling scenario as e.g. in B. 

Fig. n^ shows the experimental setup which we have 
used. The probe carrier was mounted on a precisely 
balanced horizontal linear bearing (left in Fig. nl). A 
second bearing (middle) was driven by a stepping mo- 
tor via a crankshaft with adjustable eccentricity. Both 
bearings were connected by a piezoelectric force sensor 



which allows to measure the driving force with 2 kHz 
time resolution. The motor was computer controlled, 
i.e. in precisely 2,000 impulses the motor axis revolves 
once. This high angular resolution provides quasi steady 
motion. We checked, that the finite step size per com- 
puter signal does not infiuence the convective behavior. 
The entire mechanical device was fixed on an oscillation 
damping table. The container of size 6 cm x 10 cm, con- 
sists of transparent plastic material. It was illuminated 
from the top and was monitored by a videocamera, which 
was located in a distance of 4 to in a direction perpen- 
dicular to the axis of oscillation. With this camera we 
observed the dynamical processes in the container. For 
that purpose the camera was levelled exactly at the ver- 
tical height of the granular material at rest. 
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FIG. 1. Top: Experimental setup (see text for expla- 
nation). Bottom: Snapshots of the surface region of the 
material in the container. The dashed line leads the eye 
to a fixed height above the container bottom. 

The lower part of Fig. H displays snapshots of the sur- 
face region of width 7 cm. The container is horizon- 
tally shaken with amplitude A = 0.15 cm (the same 
used in the investigations described below) and frequency 
f — 25 sec~^ . As discussed in H , we observe a small heap 
in the center. The heap-height fluctuates between about 
2.3 and 3.7 mm with respect to the level at rest. The 
dashed line is drawn 3 mm above the level at rest. The 



time dependence of the height of the surface on forcing 
parameters is subject of our interest. As shown below, 
the height fluctuates periodically with a period of about 
2.2 sec, which is 55 times the period of the driving oscil- 
lation. 

The swelling oscillation of the material has been stud- 
ied qualitatively for various combinations of amplitude 
A, frequency / and filling height and for a number of dif- 
ferent materials and grain sizes (20/im < r < 100/im). To 
analyze the effect we recorded time series of the height of 
the material surface in the middle of the container. Fig. g 
displays the height over a period of 40 sec for different 
driving frequencies. 
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FIG. 2. Height of the granular material in the mid- 
dle of the container over time for driving frequencies / 
between 17 sec~^ and 32 sec~^. 

The height of the displayed regions is 0.4 cm each. The 
material filling height (at rest) was 2.5 cm, and the av- 
erage grain size was 100 iim. For all frequencies in Fig. 
one sees that the material height varies with time. For 
small driving frequency (/ < 22 sec"^) and large driv- 
ing frequency (/ > 30 sec~^) the height changes irregu- 
larly with time, whereas for 23 sec~^ ^ / ^ 29 sec~^ we 
observe regular, almost periodic oscillation of the ma- 
terial height. This behavior can be characterized by 
the Fourier transforms of the full series of total lengths 
120 sec (Fig. 0). For / < 23 sec^^ no characteristic fre- 



quency can be observed, i.e. the Fourier spectrum con- 
tains many frequencies. 
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FIG. 3. Fourier transforms of the oscillation of the ma- 
terial height for different driving frequencies correspond- 
ing to Fig. 2. 

At / = 23 sec~^ there seems to be a sharp transition 
into another regime, where we find a characteristic fre- 
quency of the height oscillation as indicated by a peak in 
the Fourier spectrum. For high frequency / ^ 29 sec~^ 
the amplitude of swelling becomes very small and the ef- 
fect vanishes. Hysteresis of the transition point has not 
been observed. 

As seen from Fig. 0, in the region of periodic swelling 
(23 sec^^ ^ ./ ^ 29 sec^^) the swelling frequency F is a 
function of the driving frequency /. 



0.5 

7 0.4 
o 
a> 

^ 0.3 

LL 

0.2 
0.1 



X K 15mm 

D n 20mm 
A - - A 25mm 
H 1- 30mm 




15 



20 



25 

f in sec 



30 



35 



FIG. 4. Swelling frequency F as a function of the fre- 
quency of forcing for various filling heights. 



Fig. H shows the characteristic swelhng frequency F 
over the driving frequency / for different fiUing heights. 
The data shown in Fig. |j can be reproduced with good 
accuracy, ahhough it seems not to be possible to collapse 
the data using dimensionless numbers for many combi- 
nations of the system parameters. 

In the shaken container one observes convection: the 
material flows downwards close to the walls and upwards 
in the center of the container. In the interval of driving 
parameters considered here, there is an intensive convec- 
tive material flow. Viewing the container from top one 
observes that the oscillation of material height (Figs. |2|^) 
corresponds to a varying particle velocity at the upper 
surface. When the material is swelt there is an intensive 
flow, whereas the flow in the center of the surface comes 
almost to rest when the material is collapsed. A simple 
way to qualitatively visualize this effect is to take a se- 
ries of snapshots of the container viewing from top using 
a digital camera. 

In this way one produces a series of difference pictures 
by subtracting the gray scale values of the pixels of con- 
secutive snapshots. The average value AG of the gray 
scale of the pixels of the difference pictures provides a 
measure for the material flow at the surface. Fig. || shows 
the Fourier transform of AG for a driving frequency of 
/ = 25 sec~^. We find that the variation of the surface 
flow has the same characteristic frequency as the height 
oscillation (c.f. Fig. ph. For this reason we believe that 
convective motion in the horizontally shaken container 
and swelling are closely related effects. 
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FIG. 5. Material flow at the surface (central area of 
size 6.5 cm x 2 cm) over time for f — 25 sec~^. 

From the force data F(t) measured by the sensor we 
calculated the dissipated driving energy during the nth 
shaking period: 



En = 27r/A / F{t) cos (27r/t) dt . 

n/f 



(1) 



En is the energy dissipated during the nth period by the 
entire system, i.e. by the linear bearing and by the dis- 
sipative interaction of the grains in the container. The 
time series En can be analyzed by Fourier analysis sim- 
ilar to the surface flow intensity and the height of the 



material. Figure ^ shows the Fourier transform of the 
dissipated energy, again for amplitude A = 0.15 cm and 
a range of driving frequencies /. 
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FIG. 6. Fourier spectra of the dissipated energy E for 
a range of frequencies. 

Comparing Figs. O and we notice, that both figures 
agree well, i.e. the characteristic frequencies of height 
oscillation and of energy dissipation coincide. In con- 
junction with the other experimental observations this 
leads to the assumption, that the swelling-effect is essen- 
tially connected with a time-dependent convection, thus 
a time-dependent material flow and energy dissipation 
inside the granular material. 

The influence of air has to be discussed. One could 
imagine that the convective motion of sand "pumps" air 
into the bulk of the material which leads to swelling. 
When the amount of air reaches a certain extend a "bub- 
ble" escapes and the material collapses. This mechanism 
would explain the swelling effect as well as the saw tooth 
shape of the height over time curves (Fig. ||). But re- 
ducing the pressure to p = 50 Pa, which corresponds to 
a mean free path of Xair ~ 130 /im, the swelling effect 
does not disappear. Since Xair ~ 130 /im is smaller than 
the typical free path of sand grains we can exclude air as 
playing the major role in recurrent swelling. 

The origin of the new effect is not completely clear yet, 



we suggest a tentative explanation based on the measure- 
ment of dissipated mechanical energy and surface flow: 
In the horizontally shaken container one observes convec- 
tion Q, i.e. motion of the particles with respect to each 
other. To allow for macroscopic motion in a granular 
material the material has to be diluted below a certain 
density pR before (Reynolds dilatancy [|j). When the 
local density of the material is below pR the material can 
start to flow. The convection in the container implies 
shear flow, which causes further dilution and fluidization 
(e.g. iM). When the shaken material swells according 
to Reynolds dilatancy, the material becomes diluted and 
becomes at the same time less stable mechanically. At 
a certain moment the material becomes diluted to an 
extend that it loses mechanical stability and collapses. 
Then the material starts again to swell. 

Therefore, we believe that there are two competing ef- 
fects. First, due to shear the material tends to dilute. 
Second, the material has to remain mechanically stable, 
i.e. the grains on top have to be supported by the grains 
below them and these have to be stabilized by the grains 
of the next layer etc. 

Our explanation is supported by the measurements of 
the dissipation of mechanical energy in the system and 
of the surface material flow. If the material is collapsed, 
i.e. if a large part of the granular material has a den- 
sity above pR, the grains in these regions cannot move 
with respect to each other due to the Reynolds dilatancy 
effect. Hence, in these regions there is no shear motion 
and, therefore, none or just low dissipation of mechani- 
cal energy. From this consideration we conclude, that the 
collapsed material should dissipate less energy per time 
than the swelled material. This behavior coincides with 
what we found in the measurement of the dissipated en- 
ergy. Moreover our results show, that the intensity of the 
surface flow oscillates with the same swelling frequency 
F . Therefore, the experimental results on energy dissi- 
pation and surface flow support our hypothesis on the 
origin of the swelling effect. 

Other interesting questions concern the dependence of 
the swelling effect on container shape and size as well as 
on the grain material and the amplitude of shaking. We 
did experiments with several types of containers, with dif- 
ferent grain material and with various amplitudes. While 
the properties of swelling depend strongly on the details 
of the experiment we observed the effect for a wide range 
of parameters. A more detailed description of the proper- 
ties of the swelling effect will be subject of a forthcoming 



paper. 
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